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The results of classical trajectory calculations are presented for the reaction 
CIO ->■ 0 ^ Cl + 0i. This reaction is an imp'^rt.^nt step in the chlorine- 
catalyzed destruction of ozone which Is thought to occur in the stratosphere. 
Rate Citstants have been calculated for temperatures between 220 and 1000 "K. 
The calculated rate constant Is 4.3h 10”^^ exp(-191/Tl cm^ molecule ' s ^ 

and Its value at 300 °K Is 2.3 i 0.2 ' 10^' cm^ molecule"* s“ * , about a factor 
of 2 lower than the recent experimental measurements. The calculated activa- 
tion energy of 0.38 kcal/mole is In excellent agreement with the experimental 
value of 0.44 kcal/mole determined by Civne and Nip and provides support for 
their measurements. The calculations were performed using both the quasi- 
classical and phase space trajectory sampling methods. The empirical potential 
energy surface used In the calculations was constructed to fit experimental 
data for CIO, 0^ and ClOO molecules. Other Important features of this poten- 
tial surface, such as the barrier to reaction, were varied systemat ic.il 1 v and 
calculations were performed foi a range of conditions to determine the "best" 
theoretical rate constants. The present results demonstrate the utility of 
classical trajectory methods for determining activation energies and other 
kinetic data for important atmospheric reactions. 
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1 . I NTRODUCTION 

Th.. potonUal .hroat of chlor.no-cnfalv.oc1 Joplot.on of fho a. rnfosphc-rlc 
o.one layer has rocclvod conaldcraHc attc-ntlon since Molina and Rowland'- 
suRRosted that anthropoRonlc chlorof luorocarbons »ay pncvldo a slpnlflcant 
source of atonclc chlorine In the stratosphere.’ In fact, this proposal' has ^ ^ 
stimulated a tremendous research effort In aeronomv'-’ and chemical kinetics.- 
Previously, a CIX cycle had been proposed as a possible sink for stratospheric 
orone.'”-"’ but no major sources of chlorine had been Identified. Since that 
time calculations h.ave Indicated that the chlorof luorocarbons may cause serious 

depletion of the stratospheric orone laver.^'' 

The principal active chlorine specie.^ In the stratospheric CIX cvee are 
Cl .ind CIO. At hliih altitudes (40- Ml -sml , the ozone depletion Is a result of 


the cyclic process: 


Cl O 3 ► no + Om 

CIO + 0 *■ Cl + Op . 


since atomic oxvRen .and ozone are In equilibrium, the removal of either odd 
oxyRon species Is Important. At lower altitudes, the c.ntcentratlon of atomic 
oxvRon is low enouRh that the reaction (R2) Is In competition with 

CIO + NO -» Cl + NOp 

which is the dominant process below 20 km. Other CIO reactions, 

CIO + CO ♦ Cl + COp 

CIO Cl + 0 
CIO + CIO -► products 
CIO + O 3 •• products , 

are possible, but less Important at stratospheric conditions owlnR to lower 
rate constants or reactant concentrat ions. 
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Clearly the chemistry of CIO, and particularly the CIO + 0 reaction, 
plays an important role in the overall CIX cycle. An accurate determination 
of the relative amounts of CIO consumed by each of these processes requires 
knowledge of all the rate constants for reactions (R2) - (R7) at the ambient 
temperatures, which range from 210 to 270 “K for the altitudes of Interest 
(20-60 km).^’^ Unfortunately, complete reaction rate data are not available 
for these reactions. The rate constant for (R2) lias been measured in sever.al 
laboratories at 300 ®K using ultraviolet spectroscopy,'**’'^ mass spectrom- 
etry,'^*'^ and atomic resonance fluorescence detection methods.'^ The most 
widely accepted value'^ for that rate constant is (5. 3 ‘ O.R) * 10’" cm^ mole- 
cule“'s~' with the earlier workers'**"'^ giving values between 1 and 2 x 1C"". 
Two very recent measurements' have confirmed the work of Bemand, Cl^'ne, 
and Watson'^ yielding rate constants of (5.2 + 1.6) x 10"" and 
(4.4 ± 0.9) X 10"" cm^ molecule"'s~' , respectively, at 300 ®K. These 
studies'®’*^ also provide the first experimental determinations of the rate 
constants between 220 and 426 ’’K. While the observed temperature dependences 
are both small, the data are not in agreement since one study' ^ yields a 
negative activation energy. In addition, Park^® has measured the rate constant 
for (R2) between 1000 and 1500 “K in shock tube experiments and reports an 
averaged result of (7 ± 1 . 5 ) x lO"" cm^ mole"'s-' at 1250 'K. Consequently, 
the variation of the experimental data between 300 and 1250 is 

extremely small. According to collision theory, the rate constant for reac- 
tions with zero activation energy should still exhibit a temperature depen- 
dence.^' ’^2 However, the observed temperature dependence'^"^® is slightly 
less than the minimum expected on these grounds. In addition, the relatively 


make It difficult for one 


l-irge uncertainties in the measured rate constants 
to extrapolate the rate constant data to lower temperatures with confidence. 
Reaction (R3) has been studied by Clyne and Watson?^ and Zahnlser and 

Kaufman'*' who report rate constants of (1.7 ± 0.2) 10'" .and 

(2.3 i 0.8) X 10-" cm3 molecule*' s“' at 300 ''K. respectively. The amount of 
ozone depletion predicted by the stratospheric models is very sensitive to the 
ratio Ko/Kj.'^ Thus, it is of extreme importance that the temperature depen- 
dence of these rate constants be accurately known. The rate constants for 
(R4), (R6), and (R7) have been measured at 300 ®K and are 3-S orders of magni- 
tude lower than the above. ‘ 3 ’ ‘ Photodissoc iat ion cross sections for CIO have 
been determined by Johnston et al.*'' between 22SO .and 2800 A and predissocia- 
tion is known to occur at longer wavelengths. The ground state absorp- 

t ion probability is too small for photodissociation to contribute to the over- 
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all CIO chemistry in the stratosphere except at very high altitudes. • 

Pred issoclat ion is more rapid but it is still not important at altitudes 
between 20 and 50 km.^*' 

The present paper describes classical traieclory calculations of CIO + 0 
collisions which have been performed to determine the reaction kinetics and 
dvnamics of (R2). Previous studies on the H, + D‘ *' and H, + 1 ^ '7 oxchange 
reactions, among others, have shown that classical trafectory methods can be 


used to compute accurate rate constants, reaction cross sections, product 
vibrational and rotational energy distributions, etc. One requirement of 
these calculations is that the potential energy surface for the particular 
re.iction system must be known. However, for most reactions of importance to 
stratospheric chemistry, this is not the case. In the present work, an 
att«'mpt is made to construct a reasonable potential energy surface for (R2) 



based on experimental spectroscopic and thermodynamic data that are suitable 
for use in classical trajectory calculations. The development of reliable 
methods for constructing potential energy surfaces will permit the extension 
of classical trajectory calculations to a wide range of important chemical 
reactions. 

Part II of this paper describes certain physical properties of the Cl-0-^0 
system based on the available experimental data. These properties are used as 
a guide In the construction of a family of empirical potential energy surfaces 
which are presented in Part III. The set of surfaces spans a reasonable range 
of topographies for CIO 2 . and should bracket the true potential surface which 
Is not known. In Part IV, we give a brief description of the classical tra- 
jectory method used in this work. Part V contains a discussion of the results 
of the reaction cross section and rate constant calculations performed for 
reaction (R2) for the temperature range 220-1000 ®K. The calculated rate con- 
stants are compared with the experimental measurements and used to describe 
the temperature dependence of this reaction. Finally, in Part VI we summarize 
these calculations and briefly discuss the implications of the present study 
on the chealstry of CIO in the stratosphere. 

II. PHYSICAL CHARACTERISTICS OF THE Cl-0-0 SYSTEM 

In this section we siunmarize those properties of the Cl-0-0 system which 
are used in the construction of the potential energy surface for reaction (R2) 
described in section III. First the reactant (CIO + 0) and product (Cl + O2) 
limits are discussed and then the properties of the stable trlatomlc species 
are presented along with the adiabatic correlations between the trlatomlc 
molecules and the atom-diatom limits. 
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A. The CIO + 0 Limit 

f rin mnloc^ule Is knovm to dissociate to 

The ground state of the CIO nwiecuie a . 13/2 

^ri2p ,nd jhe dissociation energy 

j^round state atoms (Cl P 3/2 ^ ^2^' 

(63.427 t 0.009 kcal/mole) has been measured by Coxcn and Kamsey”*> and by 
Durle and Ramsay.^ Until recently, the detailed shape of the Rround state 
potential curye was not well known since progressions had been observed only 
for y < 2 in absorption”"- and tor 5 < v < 9 in emlssion.77« However, 
accurate ab initio MCSCFtCI potential curves have been obtained tor the X^n 
and a 7„ states ot CIO by Arnold, whiting, and hanghotf.’^ Their calculated 
ground state dissociation energy is within 2 kcal/mole ot the above experimen- 
tal value. At 300 'K. only 1.81 ot the X^" CIO molecules is in excited vibra- 
tional states (this decreases to 0.2% at 220 "K). However, the rotational 
spacings33 sre small and a broad thermal rotational distribution exists, even 
at 220 -K. The peak at 220 'K occurs for J - 12. and slgnltlcant populations 
exist for J > 30 (1.77%). These data tor CIO are summarized in Table 1. 
based on Refs. 26. 33-35. Baaco and Morse”^ have determined the - 'Uf/j 

spin orbit coupling constant to be 318 cm-' (AE - 1.364 kcal/mole). resulting 
in 90.8% of all CIO being in the x2n,/2 state at 300 ’X (95.8% at 220 ”K). 

Thus, a thermal ensemble of CIO at stratospheric temperature, consists mainly 
of molecules with excited rotational levels within the ground vibrational and 

electronic state manifolds. 

At low temperatures only ground state oxygen atoms (’Pg) exist in thermal 
systems. In the stratosphere less than 0.1% of the oxygen atoms is electron- 
ically excited and these are considered to be separate chemical species in 
most treatments of stratospheric chemistry. There are three fine structure 
states for 0 (^Pg):”-” 3Pj. the ground state. 'Pf (0.453 kcal/mole above 
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thv pround stator and thi* V. state (0.648 kcal 'no!.- hiy-'u-r In tMan 

‘I’ ). At 220 ®K, 79.5?’ of a thermal sample will be in the e.rennd fine strue- 
ture ‘it ate. 

There are six possible fine structure-state combinations for 
CIO (?n) + 0 (3p) collisions and their relative thermal collision frequencies 
are given In Table II. At temperatures up to 1000 "K tlie ground state 

is by far the most probable. However, the .,nd 

^^0^ combinations also occur with appreciable probabllltv In low- 
temperature thermal collisions. 


B. The Op + Cl Limit 

In contrast to CIO, the lower electronic states of Op have been extrcmelv 
well characterized. Three of these states are energetically accessible as 
products of reaction (R2): (T^ = 22.64 kcal/mole) and ^ r J 

(T^. = 37.73 kcal/mole). The ground state reaction Is highly exothermic : ^ 


ClO(^’n) + 0(3p^) . ciK^'P^) ^ Opf^P^v' < 13), ah; = -S4.76 kcal'mole . (R2a^ 

I p qp 


as arc the reactions which form the first two excited states of 0 


ClO(?n) + 0(3Pg) -V Cl(?P^l + 
cio(^’n) + 0(3p^) . Cl(2p^,) + o^Ot^v' 


< 7). AH“ 

’ ri0 


< ah; 

‘208 


- 3.'’ . 26 kca ’ 1 e . 

-17.2Q kcal/mole . 


(R2b) 


(K2cl 


Hov,fv..r. thore l» presei.tly no positive exper Imentnl evijenoe lor slnxlel 
oxygen formation (l.e., reaction (R2b) or (R2c)). 

l.lp»ro*, Norrish. and Thmsh” and Banco and coworkernl^‘0.‘l l,.,ve 
obnervod vibrntionally excited 0, produced by reaction (Hda). Tbev found that 
the V - 2 to v' - 7 levels (primed quantum numbers refer to product mole- 
cules and unprlmed quantum numbers refer to reactants) were equal 1v populated. 
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the ground state) and the state (0.648 kcal/mole higher In energy than 

^P^). At 220 ®K, 79.5% of a thermal sample will be In the ground fine struc- 
ture state. 

There are six possible fine structure-state combinations for 
CIO (^n) + 0 0?) collisions and their relative thermal collision frequencies 
are given in Table II. At temperatures up to 1000 "K the ground state 

far the most probable. However, the 

combinations also occur with appreciable probability in low- 
temperature thermal collisions. 


B. The Og + Cl Limi t 

In contrast to CIO, the lower electronic states of 0^ have been extremely 
well characterized. Three of these states are energetically accessible as 
products of reaction (R2): 3 e”, ^Ag (T^ - 22.64 kcal/mole) and 


(Tg - 37.73 kcal/raole).3‘* The ground state reaction is highly exothermic : 33 
C10(2n) + 0(3p ) Cl(2Py) + 02(3 e-v’ < 13). AH! 


£ = -54.76 kcal/mole , (R2a) 

2 98 


as 


are the reactions which form the first two excited states of 0., 


C10(2n) + 0(3Pg) Cl(2p^^) + O.C^gV' < 7). 


2 v «gv i /;, » -32.26 kcal/mole , (R2b) 


298 


C10(2n) + 0(3pg) ... Cl(2Py) 02(1 eJv' < 4), 


AH 


^298 * kcal/mole . (R2c) 


Howevar, there la presently no positive experimental evidence for singlet 
oxygen formation (l.e,, reaction (R2b) or (R2c)). 

Lipscomb, Norrlsh, and Thrush^^ and Basco and coworkers^ 5,40,41 have 
observed vlbt.tlonally excited 0, produced by reectlon (R2«). They found thet 
the V 2 to V . 7 levels (prined qusntun numbers refer to product mole- 
culee end unprined qusntun numbers refer to reactsnte) were equally populated. 


and that the population gradually increased from the v' « 8 to v’ - 13 
levels. Small populations were observed for v* - O.IA and there was a large 
rate of formation for O 2 v' ■ 1, They found that approximately AST. of the 
heat of reaction was channeled Into O 2 vibration. Even though large uncer- 
tainties exist in the measured Oj vibrational distribution, these data demon- 
strate that the major products of reaction (R2) are vlbrationally excited 
©2 Eg and Cl P^. Based on these measurements’* ^ and the work of Polanyl and 
coworkers, 3 who have correlated the nature of the ’’energy release" on 
potential energy surfaces with product vibrational energy distributions, one 
can conclude that reaction (R2a) proceeds along a predominately "attractive 
and mixed" potential energy surface. 

The splittings between the fine-structure states of are sufficiently 

small ( 2 cm ) to be considered zero. Some of the Important physical constants 
for O 2 are listed in Table I. The ground state of the chlorine atom has 2p^ 
symmetry with a spin orbit splitting between the and 2?^^^ states of 

2.523 kcal/raole.30 There is no experimental evidence for preferential forma- 
tion of one of the fine structure states in reaction (R2). At low temperatures, 
no other electronic states of Cl can be formed in this reaction. 

C. Triatomic CIO, 

There are two distinct Isomers of CIO 2 . One (OCIO) is a chemically 
stable free radical at room temperature whose ground state has Cjy symmetry. 

The other, e peroxy-ltke redid (clOO) of C, eywetry. Is highly reectlve end 
hee been observed only at low tenperatures . CIO + 0 colllalons must pass 

through Interswdlate geonetrlea resenbllng the structure ClOO In order to form 
Cl + ©2 products. 
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The existence of an unstable ClOO species was first proposed by Bensort 
and Buss'*** in order to account for certain observations of chlorine chain 
reactions. They estimated that ClOO has a lower total enerRy than OCIO and 
that its high reactivity Is due to a very weak Cl-0 bond. Confirmation of 
their hypothesis has come from matrix IR**^ and ESR'*- experiments. In the 
former. Arkell and Schwager have determined the equilibrium Roometrv and force 
constants. They found ClOO to have a bond angle of 110\ a weak and elongated 
Cl-0 bond, and an 0-0 bond similar to that found In 0,. Restricted Hartree- 
Fock calculations for ClOO have identified the symmetries of the two lowest 
electronic states as ^A” and (T^ < 1 Johnston and coworkers . ? 

Wagman and Garvin.***' and. most recentlv, Clyne. McKenney. and WatsoiV* have 
revised the values for the heat of formation of CIO and ClOO. The most recent 
estimates of the Cl-0 and 0-0 bond energies are 7.7 and B2.7 kcal/nh>le, 
respectively, and ClOO Is thermodynamically more stable than OCIO by 
3.1 kcal/mole. Some of the physical constants for the electronic ground 
states of ClOO and OCIO are listed in Table 111. Figure 1 is a sohemat i- 
energy diagram of the low-lying states of C10_,. 


0 . abatlc Electronic State C orrelation s 

In order to connect the at ota- diatomic and trlatomic regions of the t .O^ 
potential surface, one must consider the adiabatic electronic state correla- 
tions under the appropriate synanetry restrictions. Shuler*''^ outlined the 
basic method for determining the synmnetr ie.s of atomic and diatomic and 

D..,) fragments under lower point groups (e.g., or Cjy nonlinear triatomic 
molecules). For the present case, 6 doublet and 6 quartet states arise from 
the combination of a atom with a ’r, diatomic. The doublet surfaces are 
the only ones of importance in the present work, since the quartet surfaces 
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ct.rrelate with high energy states of CIOO.'*® Of the 6 possible states In C 
synnetry. 3 are A' and 3 are A". The first ^V' state correlate- with the Cion 
ground electronic state and the first ^A' state correlates with the low lying 
esclted state of ClOO discussed above. The energies of the 4 other states of 
ClOO that dissociate to ground state CIO + 0 are not known. Slnllarly, 

0,(3f-) + C1(^P„) correlates to ’a' + 2 ?A" (and 3 loartets). The resulting 
correlation between reactants. Intermediates and products for (R2) Is shown In 
Fig. 2. Only the correlate to ground state reac- 

tants in Cj symmetry. 5- Since the 1 2 a' and 2 2 a" states of ClOO correlate 

to the ground state of Cl + Oj (only -8 kcal/mole a'aove 1 2g"), they arc most 
likely not bound. 


Fine structure effects can also be estimated by the use of extended point 
groups51.52 to combine the orbital and spin symmetries. All polyatomic states 
and A ) become symmetry and a one-to-one correspondence between 

reactant and product fine structure states can be made by Invoking the non- 
crossing rule*3 (see Fig. The + ’P,). * lp_,, 

* ";/2 reactant combinations are most probable In thermal systems. The 


• Irst two correlate with * Vp products through the 1 2 a" and I 2g< 

states of ClOO and the third correlates with products. Of course 

these correlations hold only within the adiabatic representation. Nonadlabatlc 
effects, such as surface hopping or electronic transitions, could lead to the 


fomatlon of 02 .^** 


These thermodynamic and spectroscopic data combined with the electronic 
correlations provide a good qualitative picture of the cl-0-0 system. Much of 
this infortmitlon 1. u«td In the construction of . set of empirical potential 
energy surface, which 1. described In the next section. 


111. THE POTENTIAL ENERGY SURFACE 


Most calculations of colllsional phenomena require as Input the potential 
energy of the system for a large volume of configuration space. Since accurate 
a b initio calculations of the potential energy are not usually available, 
empirical and semiemplrical valence-bond methods (such as have been 

devised to represent potential surfaces. In the present case, an accurate 
ab initio potential surface is unavailable and the LKPS formulation is not 
suitable, since LEPS potential surfaces always favor collinear triacomic 
geometries. Bnpirlcal angular dependencies have been artificially introduced 
Into the LEPS framework by either supplying an additional angle dependent 
term^® or by modifying the standard repulsive potentials for the diatomic 
fragments. However, these approaches do not yield satisfactory results for 
CIOO. In this section an empirical model for the ClOO potential energy sur- 
face is described based on the physical properties of the constituent species 
presented in Section II. The underlying rationale employed here should be 
generally applicable to the construction of potential surfaces for other 
chemical systems for which only limited experimental and/or theoretical data 
exist. 

The potential energy surface for the ground state of CIO, has a deep well 
corresponding to the CIOO equilibrium geometry and troughs corresponding to 
separated atom-diatom geometries.^® For reaction (R2) to occur, the oxygen 
atom must impinge on the oxygen end of the CIO molecule. The reaction coor- 
dinate for process (R2a) (l.e., the minimum energy path connecting separated 
CIO + 0 and Og + Cl geometries) must traverse the stable CIOO region of con- 
figuration space and the favored direction of approach for the colliding 
reactant species will be highly nonlinear (the Cl-0-0 angle would be 


approxlnately 110*). For the present study » the ground electronic state 
potential surface has bem divided Into 5 regions, as shown in Fig. 3: two 
separated atoB-diatosi regions for CIO ’f 0 and Cl > 0 ^* labeled I and V, 
respectively; a trlatomlc region for ClOO (III) and transition regions between 
the atoB-dlaton and trlatonlc configurations (II and IV). Mathematical 
expressions are given for the potential energy in regions I, III, and V and 
an Interpolation schmse is presented to ccaabine these formulae in the transl*- 
tlon regions. 

In the separated atcmi-dlatom regions (I and V), the only contribution to 
the potential energy arises from the dlatcnaic potential curve. For computa- 
tional convenience, Morse functions (Vj^) are used to represent the ^ 113^2 ^10 
and O 2 diatomic molecules. The potential curves are constructed from the 
experimental well depth (Dg), equilibrium bond length (Rg) and vibrational 
constant (Wg) according to:^^ 

V„{R) - Dgjl - , (III-l) 


where the parameter 0 Is given by 



(III-2) 


h, c, and p are Planck's constant, the speed of light and the diatomic mole- 
cule reduced mass, respectively. Morse potentials are generally accurate near 
Rg but they can deviate considerably from the exact potential curves at other 
Internuclear separations. The calculations of Arnold, Whiting, and Langhoff^** 
Indicate that the true CIO potential curve rises more steeply at large R. 
However, only very low vibrational levels of CIO are populated under the 
conditions of the present st'idy, and the fitted potential curve Is expected to 
be adequate. On the other hand, the Morse function for O 2 exhibits a potential 


In order to prnperlv 


well that is considerably narrower than the RKR curve. 
describe the reaction dynamics, the empirical potential surface should accu- 
rately reproduce the 0-0 Interaction within the range of vibrational levels 
that might be populated by reaction (R2a) (l.OO A < R < 1.62 A for v' < 14). 

A better fit (within 0.04 A) to the RKR curve in this range is obtained bv 
setting 2 equal to 70% of the value given in Eq. (111-2). 

The triatomlc potential region III is representt'd by Mc'cse functii'ns for 
the Cl-0 and 0-0 separations and respectively) and a double minimum 

harmonic potential for the Cl-0-0 bond angle (ct). The Morse parameters are 
determined from the experimental force constants (f)"*^ according to^^ 

2 = {f/2Dp}l/? . (111-3) 

For a nonlinear trlatomic molecule the bending potential has a maximum at 
ct = 180° (the so-called inversion barrier) and two minima at the equilibrium 
bond angles and (2TT-ap). One expression for such a bending potential is 

C' + X- ” ’ (1II-4) 

where X = a - tt. This function has a maximum at X = 0 and minima at 

X = » i{-C + . OlT-5) 

The four parameters in Kq. (I1I-4) can be determined from tlu’ bending 

force constant fp, the height of the inversion barrier E*. and the require- 
ment that Y(Xq) = 0. The derivation is given in the Appendix. 

In summary, the potential energy in each of the above three regions is 
given by: 




''ui ■ ^ + V 


REFl 


(lll-6a) 

(lll-6b) 


n 


a:ul 




\ ^RF.F2 ’ 

.here V„ and V„, sy™boli^e Morse potentials parameterized for the diatomic 

and triatomic systems, respectively. «nd constants 

.e,uued to ,.ve the three potentials the cotton reierence point as shown ,n 

nn -It R 1. For this case 
Fig. 1 (separated reactants wi e 

^ ^ t.«a V * -55.83 kcal/mole. 

V = -63.59 kcal/mole and 

"""""'The transition regions are defined by 2.2 A < < ^.0 A for 

3.1 A . R,, ^ < 3.7 A ior V„. where R„_,, la the distance between Cl 

and the center-of-nass of the 0-0 fr.rpnent. The potential energy tn these 
regions is constructed by interpolating between the triatomic and atom-diatom 

, d 4 ,.«n^ro^^Pd bv a "hyperbolic t.angent-l ike 

potentials. The interpolation is controlle 

function": 

,,, , Xix - o) - ^ i lx-p|s \. (in-7) 

“(x> ■ x' + (X - p)- 2 

which switches between the two forms of the potential enerpv. x is R„„ in 
repion if and R,,.„o in region IV. P and X are parameters which fix the 
location of the center of the interpolation region and its width, respectlvelv. 
W varies between 0 and 1 and dW/dx = 0 at x - P 1 X- >7., lie this form of 
interpolation Is contlnnous. it can introduce spurious local extrema into the 
potential energy surface. To prevent this, .an extra term. 

''PERT-V'’'”' 

has been added in each transition region to smooth the potential surface. The 
use of V„.„ also permits the placement of adjustable barriers in the 

entrance and exit channels of the potential surface. For V,^. 

V - 0.5 kcal/mole for the rero barrier case when n - a, (the only one usid) 

P 


u 


and for V , V ■ ^.0, 2.'i, and 2.0 kcal/mi>lo for P„.uo ■ 0 (no barrlor^, 

I I P iv\RK 

0.5 and l.O kcal/mole, rospoct Ively , when a ■ a^,. 

Tho entire Interpolation scheme Is shown in Fig. 1 and summarized heU'w. 


The potential energy function Is first partitioned into two parts 


* 'ill ’ *^Cl-00 ^ • 


( I tl-9a) 




'b - • «ci-cx)" 


(lll-9c) 


Then V and V_ are connected bv a similar interpolat ion to complete the con- 
A n 

striiction of the potential surface: 


■ 'b • Vo ‘ = 


(Ul-lOa'l 


V - + (I - wcR^^^nVj, 4. . :.:a . r^^^^ . i.a\ ; (iii-iob^ 


V ■ V R '' 1 o\ 

A • *'00 * 


(1 11-lOc) 


Equations (Tll-lOb) and (III-‘>b) represent the potential in the transition 
re>;ions V^j and respectively. Contour plots of the resulting potential 

surfaces for a ■ ^0, 110, and 130“ are shown in I'i^. 4. The parameters used 
in the calculation of the potential ener>;y are given in T.ible IV. 

The empirical potential energy funct Ion described above should provide a 
good description of the actual potential surface in regions I, 111, and V. A 
qualitative assessment of this potential can be made in terms of the nature of 
the energy release by using the reaction coordin.ite criterion of Polanvl and 
coworkers.*** The form.atlon of ClOO from CIO + 0 involves only mixed energy 
release since R^^ decreases and Increases as CIO and 0 approach. 

I.ikewise, the formation of ClOO fri'm Cl + 0^ Involves attractive energy release 
ns the 0-0 bond length remains nearly constant during this process. Overall, 
the potential surface for reaction (R21 is classified as "mixed" and one 
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would expect a high degree of vibrational excitation In the products,'*^ ^ In 
agreement with experiment.**^ 

While this potential energy surface seems to be very reasonable for reac- 
tion (R2a), the calculated reactlor cross sections and rate constants will be 
most sensitive to the topography ot the potential surface In transition 
region II which is largely unknown. Topographical features like the height of 
the energy barrier, the 0-0 separation at the onset of ClO-0 attraction in 
CIO + 0 collisions and the variation of the potential with a near the energy 
barrier will control the fraction of collisions which reach geometries corre- 
sponding to transient ClOO complexes and continue to products. Reasonable 
estimates based on experimental observations or a range of values have been 
used for these quantities in constructing the potential surfaces. It is hoped 
that future ab initio calculations of the potential energies will identify 
which of these choices is most reasonable. The magnitude of the energy barrier 
has been varied between 0 and 1.0 kcal/mole through the choice of Vp. The 
onset of the 0-0 attraction Is governed in part by the shape of the 0-0 Morse 
potential in and by the location of the transition region II. The 

potential surface parameters given in Table IV result in an 0-0 attraction of 
1 kcal/mole at » 2.5 A for the zero-barrier surface. Great uncertainty 

lies in the selection of the rate of attenuation of the bending potential 
Y(a) as R^q and are Increased. In the above formulation, the variation 

of the potential energy with a in region II is given by 

V(j - [1 - W(Ro^)]Y(a) . (III-ll) 

The magnitude of at a particular value of R^ c.an be varied through the 

choice of boundaries for that transition region or through the choice f_ and 
E . It must be noted that the latter affects the accuracy of the potential 
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V/XUUli.NAlj rAUC 10 

OP POOR QUALTIY 


1m rci^Lo’' fll and tould ros\;lt in <*tron<ii'j 


pinduri. rid tnirp' 


nistrlSuC ions. 


!A’. r.L\SSI';V, T’sAJF.CTORY CALCULATIONS 


The trajectory of a reactive collision corrob; or to a u'lirjui.’ pith in 
phase space, proceeding from some reactant configuration t > sene product con- 
figuration. For electronically adiabatic collisiors tais path c u. uruali' be 
d-'tormined from the laws of classical mechanics.^" P'- vided that the nc'ential 
energy surface is known. The classical motion of a ' f, atom's is 

described bv 6^4^ coupled first-order differential eqva' ims. 


yff 


dP 


-i = - ^ . = 12 

dt . ’ ’ ’ 




(’•V-lal 


rn'-ib) 


The Qj and Pj are the 2^ coordinates and conjugate moner.ra, respectively, 
reeded to describe the system and K is the classioel M’-ni' t-um of 
kinetic and potential energies). This Hamiltonian, relative tj ch'’ center-nf- 
mass, is given by: 


’"2^ S’”! 


6 

V 


BC jTi p 


-*■ V(q, . q. Qr' . 


fTV-2) 


J >r a trlatonlc system, A + BC. The coordinates (Qi , , Q-jl are thi 

Cartesian components of the vector connecting atoms H n -'J C. and I'O^^ . , .i ) 

are the components of the vector between atom A and th* 'ci ter-of-mas." o*” BC. 

p„„ and u. are the reduced masses cf atom? B and C and of atom A witii 

ilC A" dC 

diatomic BC, respectively. For ^ i 3, mimerlrul integration method-' must be 
used to compute the trajectories. 


1 


/ 


In the present work a modified version of Muckerman’s quaslclasslcal 
trajectory code^l ^^ed. This computer program contains an extremely accu- 

rate numerical Integrator (11th order predictor - 11th order corrector) which 
is required to prevent the accumulation of truncation errors^'^ when long-lived 
trajectories are computed. Individual classical trajectories are subject to 
the constraints of conservation of total energy and angular momentum 

L which are used to provide Internal checks on the accuracy of the numerical 
integration. In the present work the changes in and L were typically 

less than one part in 10^ from the start to the finish of a trajectory. In 
addition, back integration of selected trajectories to return the initial 
coordinates constituted another check of the computational procedure. 

Microanonical or canonical ensembles of colliding atoms and molecules can 
be simulated by calculating large numbers of trajectories with initial values 
for Qj and Pj randomly selected from suitable distributions. There are 


several well-defined methods for selecting these initial conditions. Two of 
these, denoted the quaslclasslcal trajectory (QCT) and phase space trajectory 
(PST) sampling methods, have been used in the present work and are described 

briefly below. 


A. Quaslclasslcal Trajectory Sampling Method 

In the standard QCT sampling method the initial conditions are chosen to 
correspond to separated reactants with some relative momentum towards each 
other. The reactants are initially given classical vibration and rotation 
energies that correspond to quantum states of the isolated molecules. As 
Eqs. (IV-1) are numerically integrated, the species will first approach each 
other and then separate as products or the original reactants. The 
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original Pa(!b k, 

i'OOH QUAUIY 


quasiclassical trajectory sampling method (QCT) has been disrus'^ed In dotnll 
by Karplus, Porter, and Sharma^^ and in a rerent review by Porter. 

The initial values of the coordinates Q“ for the atom A and diatomic 
molecule BC in the (v,J) state are Illustrated in Fig. 5. The three atoms lie 
in the X-Y plane of an arbitrary space-fixed coordinate system with 
equal to one of the vibrational turning points of the (v,J) state and the 
angle of orientation, C. of randomly selected from a uniform distribution. 
Atom A is located a distance Pq from the center of mass of BC and n dis- 
tance b (the Impact parameter) from the x-axis. is a constant chosen 

such that the atom-diatom interaction potential is essentially zero at this 
separation. In the present work, p^ = 4.5 A, which results in as>Tnptotic 
ClO+0 and Cl+Oj interaction potentials of less than 0.01 kcal/mole. The square 
of the impact parameter is selected randomlv from a uniform distribution with 
a cutoff at b < bj„, where b„, is chosen such that the probability of reac- 
tion is less than 0.01 for trajectories with b > b^^. The standard error 
associated with this statistical sampling method is minimized if b^, is set 
at its minimum possible value (b^ = 2.5 A in the present work). Qu is 
incremented by 6Q“ , the product of a random fraction of the vibr.Tricnal 
period t(v) and the relative collision velocity. Thl' ensures that the 
vibrational phase will be randomly chosen even though all trajectories are 
started at the vibrational turning points. The P| are chosen so that the 
rotational angular momentum vector of BC has random orientation and magnitude 
equal to J(J + corresponding to the Jth quantum level. Tne rela- 
tive momentum between A and BC has magnitude , where 

is the collision energy and is directed along the x-axis (Pg ® Pp * 0). 
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Thus, the Initial conditions are specified for fixed v, J, numeri- 

cal integration Is terminated when one of the Interatomic distances is greater 


than Pq *^BC‘ 

The total reaction cross section, S^., for fixed v, J, and ^ Is 


defined by 


^r^^COLL’ 


J) 


2 ^'r ^^COLL’ 
^'^^COLL’ 


(IV- 3) 


where N is the total number of trajectories and N^. is the number of reac- 
tive trajectories. The standard error (one standard deviation) in this 
statistical sampling scheme is 


c - Sf(EpQ,, , V, J) 


N - N 
‘ NN ,- 


.V/2 


(IV-4) 


Typically, < O.IN, so that at least 250 trajectories are needed to reduce 
the sampling error to less than 20%. Thermal rate coefficients K(T) are 
computed from the thermally weighted reaction cross sections: 


V , J DL V , J 

^r^^COLL’ ^^^COLL ’ (IV-5) 

where j is the exact vibration-rotation partition function, j is the 
Internal energy of BC (v,J) and k is Roltr.mann's constant. The electronic 
degeneracy factor^^ has not been Included in Eq. (lV-5) and is discussed in 
Section V. The sumnuition over v and J can generally be truncated when 
Ey^ J > 4kT without introducing significant error. 

Kor most non-hydride diatomic molecules the rotational apnclnga are small 
and ijcr iTOHS sections must be compuied for many rot/it It'n/il levels, e.g., 
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for CIO rotational levels ur» to J ■ AO have approctahle populnfloti .it ^0(' K. 
In an alternate approach, a slnple set of trajectories can he det erminet' tor 
each value of v and with .1 randomlv chosen from a thermal distrl- 

hution for a given rotativ.nal temperature Using this procedure, the 

comnuted reaction cross sections are rotationally averaged .and can he defined 


as 


S/Ecoll- Tr> ■ «j' ^ ‘‘J + I)Exp(-F.j/kTp)Sr(E^p,,_. V. J) . (IV-h) 


Similarly, the rate constant is now given by 


K(T) 


A-BC 


* J ^r^^COLL’ ^R' ^COLL ^^‘COU. ’ (IV-7) 

o 

This sampling method assumes that vibration and rotation are separable, l.e., 
E\',J “ ^v ’ *'v^I’ i^hould be valid at low temperatures. 

In the present work, rotationally averaged cross sections have been com- 


puted using 500 trajectories for each set of initial conditions (E 


COLL 


. V, 


T ). For a given value of v and T , these data are then fit to a third-order 
K K 

polynomial in ^ '•^ntegrals In Eq. (TV-7) are numerlcallv ov.aluated 

using Simpson's rule. The error bounds for the rate constants are determined 
by using ^r^^COIL* ^ reevaluating K(T). 

The main advantage in using the QCT sampling method is the control one 
has in selecting the initial cor.di t ions . The computed cross sections can be 
weighted by suitable distribution functions to simulate a wide range of experi- 
mental conditions, e.g., thermal system, molecular beam experiment, etc. How- 
ever, a major disadvantage of the CCT arnroach is th.it onlv a sm.ill fraction 
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of the roraputed trajectories (u'^uallv less than 0.1) are reactive and, as a 
result, a typical rate constant calculation requires ~10,000 trajectories, A 
much more oificient procedure for thermal sys*'eTns, vliich is also used in the 
present study, is discussed below. 

r.. Pha s e Space Trajectory Sam pl ing Method 

The phase space trajectory (PST) samp’ Inq method^' Is based on Keck's 
ph.iso space theory. This .approach can b<» much more efficient than the QCT 
sampling method because the calculation of many nonreactive trajectories is 
eliminated. In the PST method, trajectories are started in the interaction 
region of configuration space, where .ill atoms are in close proximity, and 
integrated both backward in time to separated reactants, and forward in time 
to the final species which result from the collision. Most nonreactlve tra- 
jectories are never examined because ihev do not reach this interaction region. 
The details of the method are given in Kefs. 65 and 6^. 

The dimensional phase .space, C, o'" a chemical reaction system contains 
distinct regions that can be Identified as reactants 1?^ and products f?p. 

Tf a - 1 )-dimensionu I surface S is constructed to completely separate 

the reactants and products regions of phase space, then the flux of trajec- 
tories l eaving the reactants region constitutes ctn upper bound to the rate 
constant K^. The actual rate const. mt <s obtained by eliminating the con- 
tribution to the flux from any nonre.jctive trajectories which return to the 
reactants region. 

Tlio upper bound rate ci>nstant Is calculated with the assumption that the 
system is in the thermoclionilcal equilibrium, i.<*,, thermal equilibrium exists 
throughout the entire phase space, .according to 


(IV-8) 


f X/V-l 3,/^“3 

■tefT) ■ Qr' “".I ‘“'J ’ 

V „>0 

where '/p Is the velocity component normal to S and y “ ]vs| /(dS/dqj ) 
relates the differential surface area dS to the position coordinates and 
momenta. Qj^ Is the reactants classical partition function per unit volume, 
and Is defined as 


I, 


"hj/- 3 


(IV-9) 


For the 3 body case A + BC it is convenient to use spherical polar 
coordinates Rg^,. Og^,. R^.g^,. <^a-BC’ ®A-BC choose the divid- 
ing surface as S = - qj, where q^ is Rg^ or present 

study qj ■ *^00^* simplified to 

1W2 


Kp(T) 


E ZhIct" ! 

"oo J 


q*2 X 


f “»Cl-00 f 


‘"'ci-oo "ci-oo'ci-oo 




/ "cio'*cio o^Pl-v^fRcio'^wl 


t TV _1 n '» 


and the integrals in Eq. (IV-12) can be evaluated numerically using a trape- 
zoidal rule algorithm. The dividing surface is positioned to intersect the 
reactants' channel of the potential energy surface at R^^ ■ 2.65 A, the 
value of q* for which Kg is a minimum at T ■ 300 "K. 

Kg is an upper bound to the rate constant because it includes (1) con- 
tributions from those nonrcactlve trajectories which cross the dividing sur- 
face, (2) extra contributions from react. • trajectories which make multliile 
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crossings of S and (3) contributions from trajectories which originate in 
the product.*' region of phase space. To correct classical trajectorie.* 

arc computed with initial conditions corresponding randomly selected points on 
the dividing surface^® and Integrated forward in time towards products and 
backward toward reactants according to the scheme given in Ref. 65. The con- 
version coefficient ^(T) Is the fraction of "successful" trajectories (i.e., 
the fraction of reactive trajectories with the elimination of any extra con- 
tributions from those trajectories that make multiple crossings of S). Based 
on this definition, i Is a multiplicative correction to and the thermal 

rate constant, K, Is given by 

K(T) = S(T)Kp(T) (IV-11) 

The determination of the conversion coefficient compensates for the assumption 
of complete equilibrium Invoked in Fq. (TV-8)^'^ by removing all contributions 
to the flux which are only present at chemical equilibrium (contribution (3) 
above ) . There are no t rajectorl es prese nt In the non e quilibrium case which 
are absent at chemical equilibri um. The set of "successful" trajectories is 
equivalent to the set of reactive trajectories which would be obtained in a 
conventional trajectory calculation, if the initial conditions are sampled 
from classical thermal distributions (e.g., no zero-point vibrational energy 
or other features of the QCT sampling method). 

The value of K(T), as determined by Fq. (IV-11). does not depend on the 
choice of S. Thus, minimization of K^(T) results in maxlmlratlon of S, and 
minimization of the computational effort required for the calculation of 
"successful" trajectories. In other words, most of the nonreactive trajec- 
tories do not cross the optlnuil dividing surface especially if S intersects 
the crest of a barrier on the reaction coordlniite. For rviny svstems c > 0.5 
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As n result, rate ron- 


^ C C C 

and values as large as 0.9 have been computed. 

Btants and thermal reactant and product distributions can be talrulat«d from 
only 500 to 1000 total trajectories. The PST method has undergone extensive 
critical examination and been shown to be equivalent to the conventional 
classical trajectory method. However, for certain systems with large 
react.ant vibrational and rotational spaclngs the QCT method should be used to 
correctly account for the effect of quantized initial conditions. The CIO + 0 
system is ideally suited for the PST method because the CIO vibrational spacing 
is fairly small, and the rotational energy distribution is nearly continuous 
under low-temperature thernuil conditions. 

V. RESULTS 

In this section the results of the trajectory calculations using both the 
QCT and PST sampling methods for reaction (R2a) are presented. The QCT data 
are for a single CIO rotational temperature of ”100 ®K. Thermal reaction rate 
constants, reaction cross sections and activation energies are report«>d. The 
results of the PST calculations provide rate constants for the temperature 
range of 220 to 1000 “K and complement the QCT data. Reagent and product 
energy distributions have also been determined for reaction (R2a) and these 
results will be reported separately. 

A. Quasiclasslcal Trajectory Calculations 

In the present work, more than 36,000 quaslcl.nsslcal trajectories have 

been computed for reaction (R2a). These calculations used the potent 1.i] 

energy surfaces described in section III with Ebadd " ^*5. and 1.0 kcal/mole. 

dAKK 

The Initial conditions were specified for values of , between 0.03 and 

i •wwL 

6.0 kcal/mole, v ■ 0,1 and J chosen from the thermal CIO rotational 
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distribution for 300 °K. The rotationally averaged reaction cross sections, 
rate constants, and average reactant and product energies and scattering angle 
arc summarized In Table V. The variation of bond lengths and bond angle with 
time for two typical trajectories Is shown in Fig. 6. 


The rotationally averaged cross sections, ^r^'^'COLL* ~ 

determined for v = 0 at numerous collision energies between the apparent 
threshold^^ and 6 kcal/mole and for v = 1 at ~ ^ kcal/mole. 

These data are shown in Fig. 7. The translational energy thresholds for v = 0 


are zero for “ 0 ‘irid 0.5 kcal/mole and less than 0.1 kcal/mole for the 

1.0 kcal/mole barrier case. F.xamination of the trajectory data indicates that 


reactant rotation provides all the additional energy necessary for reaction at 


low collision energies (i.e., reactive trajectories), 

even though the vibrational zero-point energy (1.223 kcal/mole) could promote 


the reaction. 


Undoubtedly the classical description the threshold region is not 
completely accurate owing to the neglect of quantum mechanical phenomena such 
as tunnelling. However, the most important contributions to the rate con- 
stants come from cross sections with 0.4 < E <1.5 kcal/mole, well above 
the observed threshold values. This is seen by examining plots of the inte- 
grand of Eq. (IV-7) which are equivalent to the thermal translational energy 
distributions for reactive collisions (Fig, 8). The maxima occur at 
^COLL ~ ^ kcal/mole and the magnitudes of the integrands nearer the thresholds 
are considerably smaller. Since the cross sections for this energy range are 
large (0.4 to 1.5 A^), the neglect of quantum threshold effects should not 
introduce significant errors into the calculated rate constants at tempera- 
tures between 200 and 1000 “K. The cross sections computed for C10(v = 1) 
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are somewhat higher than the corresponding cross sections for v = 0. However, 
less than 2 % of the CIO molecules are in excited vibrational levels at 300 °K 
and reactive collisions involving vibrational iy excited CIO cun make only small 
contributions to the low-temperature thermal rate constants. 

The calculated rate constants range from 1.26 to 0.A7 x 10“^^ cm^ mole”^s~^ 
^BARR varied between zero and 1 kcal/mole. The statistical sampling 
method introduces an error of about ±15%. These values are a factor of 3 to 5 
lower than the recent experimental rate constants for reaction (R2).^'^"*^ The 
temperature dependence of the rate constant has been extracted from the 300 “K 


cross section data by an Arrhenius fit of rate constants calculated at 280 and 


320 “K. These values for are only approximate since they are based on 

the assumption that the cross sections are independent of rotational tempera- 
ture. The activation energies estimated from the QCT data are greater than 
0.5 kcal/mole for all values of ^gARR’ fair agreement with the 

experimental value of 0.38 kcal/mole determined by Clyne and Nlp.^® 

The overall thermal reaction probability, P^, is given by 


f 

P = 

T 




Nr(E 


COLL 


,v,T) 




dE 


COLL 


/•“ 

I exp(- 

Jo 


(V-1) 


^COLL^'^’^^^COLL ‘^^COLL 


where Nj. and N are determined as in Eq. (IV-5). At 300 “K 9.2% of the 
CIO + 0 collisions are reactive for the zero barrier case. This Is in good 
agreement with the experimental estimate of 10%.^^ 


S' Phase Space Tralectory Calculations 

Thermal rate constants and energy distributions have been determined for 
reaction (R2a) using the phase space trajectory sampling method. In all, over 
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28 


subtle. For example, the CIO vibration frequency mlRht be more rnpidlv 
attenuated as the oxygen atom approaches. This w.mld also result In .« fjreater 
probability of reaction at all temperature.^. In the ab.sence of more detailed 
knowledge of the CIO + 0 potential energy surface it Is difficult to determine 
the magnitude of K, to better than a factor of two by theoretical methods 
since the rate constant is very sensitive to small modifications in the poten- 
tial energy function. However, these potential surface features affect the 
preexponential factor much more than the activation energy. 

The resulting temperature dependence for reaction (R2) implies that the 
rate constant for this process is about 20X lower at average stratospheric 
temperatures than at 300 »K. According to the recent NAS report on strato- 
spheric halocarbons.l3 ^hls would mean a decrease in the predicted amount 
of odd-oxygen depletion than if a temperature independent value were used for 
K,. The amounts of stratospheric CIO consumed by reactions (R2) and (R3) would 
also be slightly shifted. It seems safe to conclude that the body of data for 
the CIO + 0 reaction is sufficiently accurate to permit a reliable determina- 
tion of the rate of stratospheric ozone depletion if all the other neces.sary 
kinetic data were available. 


VII. CONCLUSIONS 

Rate constants have been computed for reaction 

CiO + 0 -♦ Cl + 0. 


(R2) 


for tho tomperaturo range 220-1000 *K. The reenltlng r.ite constente. 

K, - 4.36 « , 0 -n exp(-„:/T) 0^3 

experimental rate date for 300 "K. The cnlculated tempereture dependence for 
these rate constent. ere In excellent egreement with the oheervetlens of dyne 
and Nlp>» end provide en Independent verlflcetlon of their me.Tsnremente. Tl,e 
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reaction rate constant at stratospheric tenperatures Is approximately 20Z 
lower than the room temperature value due to the activation enerRv of 
0.38 kcal/mole. The results demonstrate that classical tra|ectory calculations 
based on empirical potential energy surfaces can be used to determine accurate 
rate constants and activation energies of important atmospheric reactions. 
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APP ENDIX; PA R AMETERIZA T ION OF THE BENT) INC POTEN TIAL 


The double minimum function given in Fq. (IIT-4), 


Y * A(a - it)2 + 


C + (a - n) 


- + D , 


(A-1) 


is parameterized as follows. Fixing the minima at the equilibrium bond angles 
Qp and (2;t - Op) implies 


and 


(ftp - tt)2 » -c /b/a , 


(A-2) 


A(“e - + c '-h (a^ - n)2 ^ . (A-3) 

since Y(Og) = 0. The maximum occurs when a ^ v and Y - E* which yields 
a third equation: 


E* " B/C * D (A.4, 

Finally the second derivative of Y at the nlnlmum, equals the bending force 

constant f : 

B 

2Bl3(ap - 71)2 - C) 
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(A-5) 


in four u„.„ou„» «. „. c, „„, „,. „,., 

below: 


A « fgEVr . 

hAf E*3(a. - T,)** 

B » 5 

P3 


fj,(ae n)A 


(A-h) 


(A- 7) 


(A-8) 


and 


D 


- tt)2 _ 15P*] 




where r is defined by 
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Physical constants for CIO and 0^ 



cio(2n)^ 

Op(’i-)” 

OjC^Ag)*’ 

02(llg)^ 

D(j,kcal/mole 

63.31‘^(63.427) 

117.97 

95.43 

80.45 

Tp.A 

1.546‘^d. 56965) 

1.20739*“ 

1.2157 

1.22685 

Wg , cm“ ^ 

859(853.8) 

1580.211 

1509.3 

1432.6874 

WgXp.cm”^ 

6.8(5. 5) 

11.99 

12.9 

13.95008 

F 

*^v-o ’ 
kcal/mole 

1.223 

2.250 

2.148 

2.038 

4E(v ■ 1 •<- 0) , 
kcal/mole 

2.417 

4.449 

— 

— 

Be.cm"^ 

0.643^(0.62345) 

1.445572 

1.4264 

1.4004 

^J-20* 

kcal/mole 

0.767^ 

1.724 

— 

— 


From Rot. 35. unless specified otherwise. Values In parentheses 
are from Coxon, Jones and Skolnlk^^a ^„^,i,ble until 

the present calculations were completed. 

^From Ref. 34, unless specified otherwise. 

*^From Ref. 26. 

^From Ref. 34. 

^From Ref. 37. 

^From Ref. 33. 
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T’ABLF II . Reliri 

n, or 




4 0(3pj 
+0('Vj) 

+o( V^,) 

^ ^ 

,) ^o( ’p^') 

,,) +oi-Vp 

+0(3p^) 

,^^‘’^1/2 > +0(*’p j 



TABLE III. Physical 

co->stants for 

ClOO and OCIO 

tat-* K. -* -* * S t iXrX ' 'V fc. A—J ^ .JC m 

ClOO^ OCIO 


1.83 

1.47^ 

’*00’^ 

1.23 

— 

a 

110“ 

Xi 

0 

GO 

E(0,0,0) ,kcal/mole 

3.18 

3.58*^ 

u) ^ , cm” ^ 

1440. 8^^ 

945.3^ 

( 1 ) 2 , cm ^ 

407‘^ 

447.4*^ 

0 ) 2 , cm” ^ 

373 ' 

1109^^ 


^Data from Ref. 45. 

^Data from: M. C. K. Pillai and R. F. Curl, 

J. Chem. Phys. 2921 (1962). 

*^Data from Ref. 33. 

*^The corresponding force constants are 
f^Q = 9.65 mdyne/A, =■ 1,29 mdyne/A, 

^OOCl = ^oci-onci ” 

f„_ , =0.07 radyne/A. 

00-Ov.x 
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TABLE IV. Potential surface parameters. 




r-« *r 7 ■3.-X.SS s % 

Dp,kcal/mole 

7 Mr.-T-'-'V » -c X 1 

0.A-1 

s »^a A'T m. * a-< • • • 

Kp.A 


CIO (Vj) 

64.553 

2.31153 

1.546 

% 

Oj (V^) 

120.217 

2.62747 

1.20739 


00 in ClOO 

63.592 

3.3132 

1.23 


CIO in ClOO 

7.765 

3.4577 

1.83 



E*,kcal/mole 

f„,mdyne A 
o 

“e 

ClOO bend (Vj^^) 

70. 

1.04 

110“ 


p,A 

A, A 

V .kcal/mole 
P 

hi 

hv 

2.6 

3.4 

0.4 

0.3 

3.0, 2.5, 2.0^ 
0.5 

Q 

The values of Vp 
respectively. 

correspond to 

^BARR " 

and 1.0 kcal/mole 
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table V. Sunwary of QCT roauUs 


N 

Nr 

'’t 

K 


F n 0 


8^)17 

811 

0.092 

1.27 + 0. 18>‘10'^ ' 


«... . . — - —•••••“ ' 

E r. “ ^ 

'^B.ARR 


12488 

7SA 

O.O'^A 

0.78^0.13x10"^' 



0.50 

0.60 

(J) 

17 

18 

(B> 

08“ 

101“ 

<v’ > 

7.5 

7.6 

<E^> 

<fX.> 

19.36(33^3'' 

22.66(38X3 

32.43(55X3 

31.2,1(53X3 

<^> 

6.72(12X3 

5.10(011 


F -1.0 
BARR 


10444 

411 

0.030 

0.47 »0. 11x10"" 

0.77 

N.D.^ 

N.n. 

N.D. 

N.D. 

N.P. 

N.D. 

Cf- - “■**•*•' 

-1 0-1 


FnerRlea are in kcal/moln. rat. constants In cm? moloculo 

a.«U.s ,n AH .Cu-. ,.ro .U».ns,.„U.ss, 

All Jata is for 300 K. 

Vt.rnln..A fro™ rat.- constants cn.cn, .tea ns,n« K<„ (,V-7. for 
T - 280 to 320 “K. This assumos that the av.ravtoil cioss 
sections are indoponflont of rotational tomporatnro. 

^Juantltios marked N.O. were not determined. 

‘‘Quantities in parentheses are percentages of the total energy. 
‘'The experimentallv determined value for is 25 kt.il/moU 

(^45*1 , as Riven In Ref. 41. 


TABLE 

VI. PST rate 

constant data. 


T/K 

E *0^ 

BARR 

SaRR “ koal/mole‘* 

Sarr * >‘cal/iiK>Ie‘' 

I. K 

s'* 10^ ^cm^molecule” ^ s“ ^ 


220 

2.007 

0.685 

0.234 

iOO 

2.828 

1 . 286 

0.585 

500 

A. 893 

3.050 

L.901 

1000 

10.116 

7.986 

6. 305 

II. i 



220 

0.747 

0.824 

0.874 

300 

0.757 

0.815 

0.850 

500 

0.759 

0.795 

0.810 

1000 

0.679 

0.706 

0.6Q2 

III. 

K'<10* ^rm^raolecule“^s“ ‘ 


220 

1.499±0.028 

0.56510.008 

0.20510.003 

300 

2.141x0.038 

1 .048*0.027 

0.497x0.016 

500 

3. 724 ‘.0.066 

2.425^0.039 

1.54010.033 

1000 

6.869‘0.149 

5.63810.115 

4. 36 310. 130 

^Datn 

based on samples of 1000 rralector les. 


^Data 

based on samples of 500 trajectories. 



TABLE VII. Variation of the angular dependence of the potential 

A b 

energy surface. ’ 


T 

E*'kcal/mole 


KgXlO^l 

cm^molecule” ^ s" 

•-X ^ 
•1 

KxlO^' 

cm^molecule” 

220 “K 

40 


2,008 

0.736 

1.47810.040 


70 

'b 

2.007 

0.747 

1.49910.028 


70 


2.835 

0. b44 

1.82610.061 


70 


4.001 

0.597 

2.38910.113 

o 

o 

o 

7X 

40 


2.830 

0.680 

1.92410.050 

(0.43)^ 


70 


2.828 

0. 757 

2. 14110.038 
(0.58)^’ 


70 


3.993 

0.578 

2.30810.088 

(0.38)J 


70 


5.633 

0.513 

2.89010,162 


(0. 31)^^ 


1000 

"K 

40 

'b 

10.139 

0.686 

6.95610.210 



70 

fB 

10.116 

0.679 

6.86910.149 



70 


14.309 

0.502 

7.18310.320 



70 


20.605 

0.527 

10.85210. 594 

‘*Data 

based 

on samples of 500 

3* ■ t v*««r « 

trajectories 

except r*' 

> « a-.mM 

- 70. 


« fp case (1000) and E*' - 70. - fp/4 case (300). 

’'barr * 

The values of fj are given in terms of f^ « 1.04 mdyne A, the 
experimental bending force constant for CIOO. 

The values in parentheses are the activation energies in kcal/mole 
as determined from the 220 and 300“ data. The lower limit (hard 
.sphere case) is 0.26 kcal/mole. 
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figure captions 


Fig. 1. 
Fig. 2. 


Fig. 3 


Fig. 


Fig. 


Energy diagram for the ClOj system. The energies are In kcal/mole 
relative to separated CIO + O(^Pg). 

Adiabatic electronic correlations between low-lying states of CIO 4- 0. 
CIOO. and Cl + Oj. Electronic states of ClOO denoted by dashed lines 
have not been observed experimentally or theoretically. 

Schematic partitioning of the ground state CIOO potential energy sur- 
face for a - Oj. Regions I and V correspond to sep.irated atom- 
diatom conflgur..tlons and region 111 corresponds to trlatomlc config- 
urations. The transition regions arc II .and IV. Is determined 

from for Rp^, > 2.2 A and V„ Is determined by an interpola- 
tion between and for R„, < 3.0 A. The potential In 

region IV Is specified by Vj for 3.1 < Rpi.po ' '' 

potential in region II la specified by an Interpolation between 

and V . The reaction coordinate for CIO + 0 Cl + Oj Is also shown. 
B 

The barrier is located at 0 and the CIOO cqulUbriutn geometry is 
located at o. 

Contour plots of the ground state CUW potential energy surf.ico for 

F *0 The energy contours are In kcal/mole relative to stpa 

BARR 

rated CIO + 0. The bond angle a is fixed as follows for each plot: 
(a) a - (b) a = a,. ^ 110». anu (cl a » 130». The CIOO enuilib- 

rlum geometry is located at 0 in (b) . 

The parameters used to specify the initial coordinates for the quasl- 
classical trajectory sampling method. is the displacement of 

atom A along the -x axis to effect randomization of the initial 
vibrational phase of BC. 
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Fig. 6. Typical reactive trajectory. (a) Variation of bond lengths and bond 
angle with time for a trajectory with - 1.0 kcal/mole, 

^BARR “ encounter time Is 4.2 x 10“13 s. (b) Path of this 

trajectory In space superimposed on the potential surface 

with a - a^. Selected points labeled a through e are shown In 
both (a) and (b). 

Fig. 7. Rotatlonally averaged reactive cross sections for CIO + 0 ^ Cl + Oj 
at 300 “K. The points denoted by 0 and o are for v = 0 and 

V » 1, respectively. The smooth curves are polynomial least square 
fits to the V = 0 data. Statistical error limits are shown for the 

V - 0 data. The results are for: (a) . 0. (b) . 0.5. 


and (c) » 1.0 kcal/mole. 

Fig. 8. Translational energy distributions for reactive trajectories. 

*"^^C0LL^ * V = 0, T - 300 “K). The poUmomlal least squares 

fits to the rotatlonally averaged cross sections were used. The 

dashed curve Is the distribution for all collisions. F(F ) 

'COLL ' * 

Fig. 9. Arrhenius plot of the rote constant data for CIO + o * C1 + 0,. The 
solid curves are calculated P!!T rate constants. The da.shed curve Is 

the "best" PST result for E_ « 0 and f = f /-) tu j 

BARR ^ ‘b points 

represent experimental measurements: (a) Bem.ind. Clyne and Watson, 

(b) Clyne and Nlp‘f (....). (e) Zahnlser and Kaufman'’ (— ), 

(d) Freeman and Phillips."' (e) B.,sco and Dogra.'’ and (f) Park.2« 

Fig. 10. Contours of equal potential energy as a function of R and a 

CIO 

**00 **1 ^ fpF ^BARR " *** contours are for V « 1, 

5 and 10 kcal/mole. The dashed Hues are for R^,^ . r^ and a . a^. 

<o) E. . 70 kcal/mole and (• - f^, (h) r. . ,0 kcal/mole and 




fB/4. 
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